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Abstract 

_j_ , Available experimental data on decay rates for B — ► ipK, and decay rates 

£> \ and transversity amplitudes for B — * ipK* are used to investigate nonfactoriza- 

tion contributions in these decays using five different theoretical models for the 
^-j- . formfactors. Using the knowledge on nonfactorization so gained, we study the 

•/^ processes B s — ► ip(r),rf), B s — ► ifxj) and B — > r ip{2S)K* . We find that present 

experimental data for the last two processes are consistent with the predictions 
■ of most of the models considered. 

&\ (PACS numbers: 13.25.Hw, 14.40.Nd) 

£ ; 1 Introduction 

j> ! It was shown in [T], 0| that factorization approximation used in conjunction with 

formfactors derived in most commonly used models failed to account for the ratio 
(B — > ipK)/(B — > ipK*) and the longitudinal polarization Tl/T in B — > ipK* decays. 
Subsequently, it was realized that nonfactorized contributions could play an im- 
portant role in these color-suppressed decays, and it was demonstrated 0, ||, |6| how 
such contributions could lead to an understanding of B — > tpK* and tj)K data. 

Our aim in this work is to investigate the role of nonfactorization and, where 
relevant, final state interactions (fsi). The most recent CLEO data enable us a 
complete amplitude analysis of B — > tpK* decay. We use this to determine the three 
partial wave amplitudes, S, P and D, and the two relative phases. We exploit this 
knowledge in our work (see || for a more restricted analysis). 

With the knowledge gained from the study of B — > ipK and ipK* decays, we have 
also investigated the processes B s — > ipr], ipr]', ipcj) with special emphasis on the role 
of nonfactorization. Furthermore, we have extended our analysis to channels with 
i/j(2S), instead of ip, in the final state. 
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2 (B, B 8 ) ip(il>(2S))P Decays 

We begin with the formulation for B decays involving ip (or ip(2S)) and a pseudoscalar 
particle in the final state. The decay amplitude, in the notation of [[5], is written as 



A(B^^P) = (^P\nf\B) 
Gp 



v*v 

V2 <* 03 



a 2 ^P\{bs){cc)\B) + d {^P\H^\B) 



where the brackets (bs) etc. represent (V — A) quark currents and 

?4 8) = lj:^ a s)(cX a c) } (2) 

^ a 

«2 = ^ + C 2 - (3) 

N c is the number of colors (taken to be 3) and A a are the Gell-Mann matrices. C\ 
and C2 are the standard Wilson coefficients for which we take the values ||, 

C x = 1.12 ±0.01, C 2 = -0.27 ±0.03, (4) 

which are consistent with the choice in J8|. 

While the second term in ([]]), the matrix element of Tiffl, is nonfactorized, the 
first term receives both factorized and nonfactorized contributions |§. We introduce 
the following definitions to proceed further, 

(V|(cc)|0) = e"m^, (5) 

(P\(bs)\B) = (pB+pp-^^q) Fr^ + ^^q^tf). (6) 
\ * / fj, Q 

The first matrix element in (Q) is then written as 

(ipp\(bs)(cc)\B) = (P|(5 s )|B><v«|(ac)|o> + <PVI(5s)(a:)|s>»/ (7) 

= 2m^UC P (e.p B ) [if V) + F^ nf (q 2 ) 
where we have defined 

(i;P\(bs)(cc)\B)^ = 2m^UC P (e.p B )F? )nf (q 2 ). (9) 
In (H) and @, Cp has the following values, 

sjl (cos 6 P + -^sm9 P ) P = r\ 
C P = \ cos e P - sin e P ) P = i (10) 

1 P = K°,K + 
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with the 77 — 77' mixing angle Op = —20°. The superscript (1) in F[ {q 2 ) 
'color-singlet'. The nonfactorized matrix element of Tiffl is parametrized as 

^P\Hf\B) n f = 2m^UC P (e.p B )Fl 8)nf (q 2 ). 

In terms of the definitions in (4) - (8), the decay amplitude for B 



denotes 



as 



where 



A(B -> ipP) 



Gp 
V2 



V^V ca al ff 2m^UC P (e.p B )Fl 



BP ( 2\ 



a 2 



a 2 



^(8)n// 2 n 



F/^g 2 ) a 2 F/^g 2 ) 



(11) 

■0-P is written 
(12) 

(13) 



As a short-hand notation, we introduce two parameters XFi and t^p 1 as measures 
of nonfactorized contributions, as follows: 



a 2 ff = a 2 ( 1 



a 2 



where, 



(8)n// 2 > 



XFi = 



+ 



(14) 



(15) 



d F? p {q 2 ) F? p {q 2 ) ' 

Note that as C\ja 2 is of the order of 10, the nonfactorized contribution from color- 
octet current is greatly enhanced. A departure of £p x from unity, or xf% from zero, 
signals nonfactorized contribution. 

In terms of the quantities defined in (H) - dl3|) , the decay rate for the exclusive 
channel B —>■ ipP is given by, 



T(B^^P) = ^\V cb \ 



Wrj 



eff 

a 2 



fi\C P \*\k\* ™ 



(16) 



where |fc| is the momentum of the decay products in B rest-frame. 

Other parameters we used were, V cs = 0.974, V cb = 0.04 0, = 0.384 ± 0.014 



GeV and /W 2 <j) = 0.282 ± 0.014 GeV |TT 



2.1 B ip(ip(2S))K decays 

We first consider the decay Z? + — > -^iV 1- whose branching ratio is more precisely 
measured than that of the neutral mode [[HJ . The decay rate formula for this process, 
(|16|), can be rearranged and written as 



/no Ar>r, _ m19 /-I _T r-2 1M/2 It/ I It/ l_ X T^BK I 2 N ' V 1 '/ 



(78.422 x 10 12 GeV" 2 sec" 1 ) 1 ^ |T/ CS | IV^/,^ 
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Table 1: Model predictions of formfactor Fi(q 2 ) at q 2 = or mL 2 n. In CDDFGN 
model, rj stands for 7] s , the octet member. This scheme cannot handle rji, the flavor 
singlet. 





BSW I 


BSW II 


CDDFGN 


AW 


ISGW 


B + - 




0.565 


0.837 


0.726 


0.542 


0.548 


B+ - 


-> ip(2S)K + 


0.707 


1.31 


0.909 


0.678 


0.760 


B° s ~ 




0.49 


0.726 


0.771 


0.534 


0.293 


B° s - 




0.613 


1.14 


0.964 


0.668 


0.475 


B° s - 




0.411 


0.609 




1.06 


0.463 


B° s ~ 


-> V(2SW 


0.514 


0.954 




1.33 


0.752 



A departure of ^ from unity signals the failure of the factorization assumption 
for a particular model-value of the formfactor. We employed five different models 
for F® K {mfy. They were: {%) BSW I |T2|], where the formfactors are calculated at 
q 2 = and extrapolated using a monopole form with the pole masses given in ||12|| , 
(ii) BSW II, where a dipole extrapolation is used for Fi(q 2 ), A 2 (q 2 ) and V(q 2 ), with 
the same pole masses as in |12| |, (Hi) CDDFGN [|13j , where the normalization of 



the formfactors are extrapolated using a monopole form, (iv) AW ||14|| , where the 



We allowed 

F BK 

(rnV) to vary continuously and determined the allowed values of 



formfactors are evaluated at the zero-recoil point corresponding to the maximum 
momentum transfer and then extrapolated down to the required q 2 using a monopole 
form, and (y) ISGW where the formfactors are calculated at the maximum q 2 
and extrapolated down to the needed value of q 2 with an exponential form. The 
predicted formfactors in these five models relevant to the processes of interest are 
shown in Tables [1] and | 

from data. The results are shown in Fig. [I], where the dots represent the values of 
F® K (m\) in various models which read from the left, AW, ISGW, BSW I, CDDFGN, 
and BSW II, in that order. ^ different from unity (or XFi different from zero), signals 
presence of nonfactorization contributions. 

We repeated the above analysis for B° — > ipK° and B + — > ip(2S)K + . The results 
are displayed in the plots of Fig. 0. For B + — > ipK + , this is simply another way to 
display the results shown in Fig. |1|. In Fig. ^ we have plotted the branching ratios 
predicted in the five different models we have considered against the parameter 
One can read-off the amount of nonfactorization needed to understand the measured 
branching ratios in each model. Clearly, nonfactorized contributions are needed to 
explain data. For example, B + — > ipK + branching ratio requires that the nonfactor- 
ization parameter £ Fl is in the range (2 - 3.5), while B° — > ipK° data require ^ to 
be in the range (1.5 - 3). BSW II model requires the least amount of nonfactorization 
due to the fact that a dipole extrapolation of the formfactors allows the factorized 
term to be larger thereby reducing the necessity for the nonfactorized contribution. 
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Table 2: Model predictions of Ax(mV), A 2 (m3,), and V{rri^) formfactors for the pro- 
cesses B -> ^.fiT*, -> ip(2S)K* and ff s -> ^0 





A, 




v 


X 


y 






BSW I 

UU V V J- 


0.458 


0.462 


0.548 


1.01 


1.19 






BSW II 

J. J U V V J- .L 


0.458 


0.645 


0.812 


1.41 


1.77 


B - 




CDDFGN 


0.279 


0.279 


0.904 


1.00 


3.24 






AW 

A V V 


0.425 


0.766 


1.19 


1.80 


2.81 






ISGW 


0.316 


0.631 


0.807 


2.00 


2.56 






BSW I 


0.549 


0.554 


0.685 


1.01 


1.25 






BSW II 


0.549 


0.924 


1.27 


1.68 


2.32 


B - 


+ ^(25)^* 


CDDFGN 


0.334 


0.334 


1.13 


1.00 


3.39 






AW 


0.509 


0.916 


1.49 


1.80 


2.94 






ISGW 


0.438 


0.875 


1.12 


2.00 


2.56 






BSW I 


0.374 


0.375 


0.466 


1.00 


1.24 






BSW II 


0.374 


0.523 


0.691 


1.40 


1.85 


B s - 




CDDFGN 


0.265 


0.279 


0.919 


1.05 


3.47 






AW 


0.449 


0.703 


1.34 


1.56 


2.98 






ISGW 


0.237 


0.396 


0.558 


1.67 


2.35 




Figure 1: Allowed region (bounded by the two curves) of £p 1 as a function of 

defined by B + —>■ K + ip. The dots show the model predictions of the formfactors; from 

left to right: AW, ISGW, BSW I, CDDFGN, BSW II. 



5 



0.0015 



0.001 



55 

03 



0.0005 




■ 0.00i5 




0.00/5 



0.007 



55 

a? 



0.0005 




BSWI 

BSW11 

CDDFGN 

AW 

1SGW 



Figure 2: Branching ratios as functions of £ Fl in each model. Horizontal lines define 



the branching ratio bounds to one standard deviation. Data from [10| for B — > ipK 



and p[ for B -> ^(25)7^ 



2.2 S s -> ^(^(25 , ))y7,77 / Decays 

The analysis of the branching ratio data for B — > tp(tfj(2S)K decays can be used to 
predict the branching ratios for B s decays into ip or ip(2S) and rj or rf if we assume that 
the amount of nonfactorized contribution is approximately independent of the light 
flavor. The calculation proceeds in a straight forward manner. We show the results 
in Fig. |3] where the branching ratios for B s — > (iprf), (i/)(2S)t]), (iprj') or (tfj(2S)i]') are 
plotted as functions of £ Fl for the five models we have considered. 

For a given £ Fl , BSW II model produces the largest branching ratios for B s — > 
if)(ip(2S))r]. For B s — > ip(ip(2S))r]' decays, the largest branching ratios for a given £ Fl 
are generated in the AW model followed by those in BSW II scheme. In order to get 
some feel for the predicted branching ratios, we have presented the model-averaged 
branching ratios in Table [5] for £ Fl = 1 (factorization), £ Fl = 2 and £ Fl = 3. 

3 (5, B s ) -> ^{2S)V Decays 

The decay amplitude for B — > ipV in the notation of || is 

A{B^W) = ^V cb V; s a 2miP U (18) 
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Figure 3: Branching ratios as a function of fa t in each model. In CDDFGN model, r\ 
stands for r/g and there is no prediction for rf . See Fig. ^| for legend. 

Table 3: Average branching ratios predicted by the theoretical models for three 
choices of fa. 





(Factorization) 
xlO^ 3 


(&i = 2) 
xl0~ 3 


(fo = 3) 

XlO- 3 


B(B S - 




0.037 


0.15 


0.33 


B(B S - 




0.017 


0.07 


0.16 


B(B S - 


-> %j)rf) 


0.084 


0.34 


0.76 


B(B S - 


- ^(2S)t/) 


0.032 


0.13 


0.28 



x [{m B + m y )(e* 1 . e * 2 ) (^(mj) + A?^ + ^ 4 8)n/ ) 

(e* 2 .(PB ~ Pv))(et.{p B + Pv)) f A BV (rn 2 , , M)nf C X 4 (8)n/ N 



(m B + m y ) w 1 V v ^ a 2 

In fllS]), ei and e 2 are the ^ and V polarization vectors respectively. Af v , A 2 V and 
V BV are the formfactors defined in [0 which contribute to the factorized part of the 
decay amplitude. A x n ^ , A^™^ and V^ n f are the nonfactorized contributions arising 
from the color-octet current products (cc)(sb), the analogues 

of pQW in (p). And 

Af )n/ , A 2 1)n/ and \/« n/ are the analogues of in ©. 
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The following definitions 
the equations that follow. 
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1 + 



1 + 



A 



(l)nf 



+ 



in analogy with Eq. (|I3D, facilitate shorter forms for 

d Aj 8)w/ 



« 2 Af V ( 



1 + 

a 2 



1,21 



y(l)n/ ^ y(8)nf 

V BV (ml) + V 2 V BV (ml 



a 2 



-Xv- 



(19) 



B — > -01/ decays can be discussed in three equivalent basis-amplitudes: Helicity 
basis (H , H + , HJ), Trasversity basis (A ,A\\,A±), and Partial-wave basis (S,P,D). 
They are related through the following definitions [IT], [T8], 



# = 






+ = 




1 


= 


4=5- ip + 
73 ^2 


i 

76 



(20) 



^0 


= H = - 


7S s + i 






= 7f<^ 


+ -»_) = 






= 7f ( ^ 


-#-) = 


P. 



(21) 



5 = 


\S\e i5s 


p = 


\P\e i5p 


D = 


\D\e iS ° 



All amplitudes in ( P0|) and (]2lD are, in principle, complex, their phases defined by 
the following (total angular momentum J = 0, uniquely determines the spin angular 
momentum once the orbital angular momentum is specified): 



(22) 



The process of generating complex amplitudes in terms of 5$, Sp and is as 
follows: Helicity amplitudes are evaluated directly from fll8|) before final-state in- 
teraction (fsi) phases are put in. This allows us to determine the real partial wave 
amplitudes, before fsi, through 
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P = ±={H+-H.) = A ± (23) 

D = A=(H + + H_+2H ) = 

Once the real S, P and D amplitudes are determined, their phase are put in by hand 
as in (]2"2]), and with these complex S, P and D wave amplitudes one can write down 
complex (H , H + , H_) from (pOj) or complex (A ,A||,Aj_) from (pT|). 

With these definitions the following expression for the decay rate are obtained 

r(s->#0 = tt^t(\ h o\ 2 + \h + \ 2 + \h_\ 2 ) 



\k\ 



B 



\k\ 



-(|A | 2 + |A||| 2 + |AJ 2 ) 
-(|S| 2 + |P| 2 + |P| 2 ), (24) 





|#0 


2 




2 + l 




2 + 




2 



(25) 



where k is the momentum in B rest-frame. 

The longitudinal and transverse polarizations are defined as, 

P Tl 

Pt = 1-Pl. (26) 

We note from fl2"U| ) that both P L and Pt depend on the relative phase between the S 
and D waves, (5$d — $s ~ $d), through cos 5sd m a compensatory manner such that 
Pl j tPt = 1 • Alternatively, one can define in term of transversity amplitudes ]17 , [F| , 





4l 


2 


\A 


2 + 


4l 


2 + 


A ± 


2 



(27) 



P± = ^ = ^— L — ^ — (2s) 





A ± 


2 




2 + 


*\\ 


2 + 




2 



r 

P = (29) 

From one notes that while Py depends on Pi is independent of the 

strong phases, 5s, 5p and 5d- Eq. ( ^0|) also shows that the only way information on 
P wave phase can be obtained is via r + oc |P+| 2 or T_ oc |P_| 2 . IV = T + + T_ is 
independent of the P-wave phase. With these definitions, we now consider the specific 
case of B — > ipK* . 



3.1 B —¥ if)K* Decays 

With the definitions introduced in the preceding section, the helicity amplitudes for 
B — ► ^fT* decays are (before fsi phases are introduced), 

G 

Ho = -—^V cb V* s f^m^(m B + m K *)a 2 Af K \m 2 b )(a£ 1 - b&x) (30) 
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where |fl| 



Gp 

H± = --j^/ ch V2J^m^{m B + m K *)a 2 Af K *(ml)(£,i =F c£ v y) (31) 



2\k\ 2 m 2 B 
m^m K *{mB + ttik*) 2 ' 

2\k\mB 
(m B + m K ,) 2 ' 
2 1 



(32) 



x 



A? K *(miy 
V BK \m\) 
A? K *(m: 



2 - 



With these helicity amplitudes, we can write the transversity amplitudes or the partial 
wave amplitudes via (|2~l"D and (f23|). 

We start with an amplitude analysis using the latest CLEO data 0]. Working in 
the transversity basis, they @ determined the following branching ratio, dimensionless 
amplitudes (denoted by A , A± and A\\) and their phases, 

B(B — > ^if*) = (1.35 ± 0.18) x 10" 3 

M 12 

|A | 2 = L_2i_ = 0.52 ±0.08 
\A I 2 

|AJ 2 = ^^i- = 0.16 ±0.09 (33) 
I4J2 

|in| 2 = L-lL = 1 - |Aj 2 - |A | 2 = 0.32 ±0.12 

011 = 3.00 ±0.37 
<j>± = -0.11 ±0.46, 

where 0m and 0j_ are the phases of the amplitudes An and A±, respectively, with the 
choice 0o = 0. The important feature is that the amplitudes are relatively real. 

Using the relation between the helicity, transversity and partial wave phases, the 
CLEO analysis (1331) can then be restated in the following equivalent forms, 



Helicity basis: 
\H I 2 

\H \ 2 = LJ1L = 0.52 ±0.08 

Irr |2 

lii+l 2 = J_±L = 0.014 ±0.034 (34) 



10 



l-|#o|-|#4 
2.92 ± 1.70 
3.01 ±0.29, 



0.47 ±0.08 



Partial-wave basis: 



I op 

|S| 2 = ^f- = 0.77±0.12 
\P\ 2 

\P\ 2 = l -±- = 0.16 ±0.09 (35) 

F 



\D\ 2 

\D\ 2 = l —±- = 1 - \S\ 2 - \P\ 2 = 0.073 ±0.044 



3.07±0.19 
P = -0.11 ±0.46 
D = 0.17 ±0.44, 

The phases in ( p4]) and ( |35"D are evaluated relative to </> as in the CLEO analysis 
From (|35|) we note that |5| > |P| > |-D| as one might intuitively anticipate. 

Next, having determined the helicity amplitudes, we are in a position to extract 
information on the parameters £i, £2 and £y in a given model for the formfactors. We 
summarize our method below. 

The helicity amplitudes in terms of £1, £2 and £y are given in Eqs. (|30"D and (|31|). 
The parameters £1, £2 and £y, representing nonfactorization, are then obtained from 
the constraints, 

l^ol 2 = ? ( ^ 2 ~? X6l)2 I , , 9 =0-52 ±0.08 (36) 

l#+| 2 = TTTi ^fc^r , , , 9 = 0.014 ±0.034 (37) 

1 +l 2(1 + cV&i) + (a - &x£ 21 ) 2 v ; 

and 

b{b^^k*) = ^IJ^|v; fe | 2 |v;//X(^ + ^) 2 «2l^iK)l 2 ei 2 

x{(a-6x£ 2 i) 2 ±2(l + cV£ 2 /1 )} 
= (1.35 ± 0.18) x 10" 3 (38) 

where we have defined the ratios, 

£21 = t 2 - and £ V i = t~- (39) 
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-2 -1 1 2 3 4 5 



Figure 4: The region in £21 and £vi plane allowed by the latest CLEO II measurments 
of \Hq\ 2 (vertical bands) and \H + \ 2 (horizontal band) for B — > ^ii"* in BSW I model. 

The regions in £21 and £yi space that explain the experimental value of polarization 
(l-ffol 2 ) are shown by the two vertical bands in Fig. |], while the region between the 
two horizontal curves corresponds to non- negative values for \H + \ 2 within the error. 
The boxes in Fig. || show the four solutions we get by solving (|36|) and ( p7|) for £ 2 i 
and The errors in £ 2 i and £vi are correlated as parts of the boxes lie ouside the 
overlap of the horizontal and vertical bands. 

Clearly, in BSW I model, within errors, there are solutions with £ 2 i = 1 and 
£yi = 1 i.e. £1 = £2 = £v- This is the class of solutions discussed in ||. The value of 
£1 is then obtained from ( |38| ) and does not allow £1 = 1 solution. We repeated this 
procedure for other models of the formfactors. the results are tabulated in Table 4. 

From Table 4 it is evident that only BSW I model permits solutions with £1 = 
£2 = £v (but £1 7^ 1); other models do not allow solutions with £1 = £2 = £v- 
There are no solutions to the latest CLEO data [[F] consistent with factorization i.e., 
£1 = £2 = = 1 in any of the models we have considered. 

We chose to work with CLEO data only rather than use the world average for 
the longitududinal polarization because it is only the CLEO data which allow a 
complete determination of the decay amplitude. For the record , the world average 
(our estimate) of all measurements [19|, ^C], |7j of the longitududinal polarization 
in B — > ipK* is 

P L = 0.66 ± 0.05. (40) 
To one stadard deviation there is no overlap of the world average and the CLEO data. 
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Table 4: Solutions of £1, £21 and £vi using the latest CLEO II measurements of 
branching ratio, polarization and \H+\ 2 for the process B — > ipK*. In the table only 
positive solutions for £1 are shown since for every solution of £1 there is another which 
is its negative. 



Model 


Solution 1 


Solution 2 


Solution 3 


Solution 4 


BSW I 


61 


1.79 ±0.62 
1.03 ±0.29 
1.36 ±0.61 


1.79 ±0.62 
3.77 ±0.29 
1.36 ±0.61 


1.26 ±0.55 
0.46 ±0.68 
2.74 ± 1.20 


1.26 ±0.55 
4.34 ±0.68 
2.74 ± 1.20 


BSW II 


61 


1.79 ±0.62 
0.74 ±0.21 
0.92 ±0.41 


1.79 ±0.62 
2.70 ±0.21 
0.92 ±0.41 


1.26 ±0.55 
0.33 ±0.49 
1.84 ±0.83 


1.26 ±0.55 
3.11 ±0.49 
1.84 ±0.83 


CDDFGN 




2.94 ± 1.00 
1.04 ±0.29 
0.50 ±0.22 


2.94 ± 1.00 
3.80 ±0.29 
0.50 ±0.22 


2.07 ±0.91 
0.46 ±0.68 
1.01 ±0.45 


2.07 ±0.91 
4.38 ±0.68 
1.01 ±0.45 


AW 




1.93 ±0.67 
0.58 ±0.16 
0.58 ±0.26 


1.93 ±0.67 
2.11 ±0.16 
0.58 ±0.26 


1.36 ±0.60 
0.26 ±0.38 
1.16 ±0.52 


1.36 ±0.60 
2.43 ±0.38 
1.16 ±0.52 


ISGW 




2.60 ±0.90 
0.52 ±0.15 
0.63 ±0.28 


2.60 ±0.90 
1.90 ±0.15 
0.63 ±0.28 


1.83 ±0.80 
0.23 ±0.34 
1.27 ±0.57 


1.83 ±0.80 
2.19 ±0.34 
1.27 ±0.57 
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Table 5: The branching ratios, \Ho\ 2 and \H-\ 2 for the process B s 
£21 and £yi the values in Table f§ calculated for B — > ^if*. 



-00 using for £1 



Model 


Solution 1 


Solution 2 


Solution 3 


Solution 4 


BSW I 




BR xlO" 3 

#o| 2 
#_| 2 


0.91 ±0.14 
0.50 ±0.07 
0.49 ±0.08 


0.75 ±0.12 
0.39 ±0.08 
0.60 ±0.09 


0.88 ±0.14 
0.48 ±0.07 
0.50 ±0.08 


0.77 ±0.12 
0.40 ±0.08 
0.58 ±0.10 


BSW II 




BR xl(T 3 

Ho\ 2 
H„\ 2 


0.91 ±0.14 
0.49 ±0.07 
0.49 ±0.08 


0.75 ±0.12 
0.39 ±0.08 
0.59 ±0.09 


0.89 ±0.14 
0.48 ±0.08 
0.51 ±0.08 


0.78 ±0.12 
0.41 ±0.08 
0.58 ±0.10 


CDDFGN 




BR xl(T 3 

H \ 2 
H_\ 2 


1.20 ±0.19 
0.47 ±0.08 
0.52 ±0.09 


1.20 ±0.19 
0.47 ±0.08 
0.52 ±0.09 


1.20 ±0.19 
0.46 ±0.08 
0.52 ±0.09 


1.20 ±0.19 
0.46 ±0.08 
0.52 ±0.09 


AW 




BR xl0~ 3 

#o| 2 
H„\ 2 


1.70 ±0.28 
0.53 ±0.06 
0.45 ±0.08 


0.96 ±0.17 
0.19 ±0.08 
0.79 ±0.08 


1.50 ±0.27 
0.49 ±0.07 
0.50 ±0.06 


1.00 ±0.18 
0.24 ±0.09 
0.73 ±0.13 


ISGW 




BR xl0~ 3 

#o| 2 
H_\ 2 


0.82 ±0.15 
0.57 ±0.06 
0.41 ±0.07 


0.42 ±0.07 
0.15 ±0.08 
0.80 ±0.08 


0.71 ±0.14 
0.54 ±0.06 
0.46 ±0.06 


0.43 ±0.07 
0.23 ±0.11 
0.76 ±0.13 


Experiment |2L 22 


B(B S -> ^0) = (0.93 ± 0.33) x 10~ 3 
P L = \H \ 2 = 0.56 ±0.21 



However, within errors there are values of £21 and £yi, shown in Table |] that fit the 
world- averaged Pl- 



3.2 B s — > ifjcf) Decays 

If it is assumed that that nonfactorization contributions are independent of the flavor 
of the light degree of freedom we can use the values of £1, £2,1 and £y t i in Table f| 
to predict the branching ratio, polarization and transversity for B s — > iptfi. Indeed, 
the branching ratio, longitudinal polarization and trasverse polarization (|if_| 2 ) for 
B s — > ijxj) can be related directly to the those for B — > ipK* by eliminating £1, £21 and 
£vi. The result is shown in Table |5|. The following experimental data for B s — > ifxp 
are now available IplL 1221 



B{B S 
Pl 



(0.93 ±0.33) x 10~ 3 
0.56 ±0.21. 



(41) 



By studying the results presented in Table [|, we note the following: (i) Most 
of the predictions (those given by BSW I model, BSW II model, CDDFGN model 
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Table 6: The branching ratios, \Ho\ 2 and \H-\ 2 for the process B 

ipK* 



ip(2S)K* usin^ 



Model 


Solution 1 


Solution 2 


Solution 3 


Solution 4 




BR xlO" 


-3 


0.80 ±0.19 


0.43 ±0.11 


0.62 ±0.17 


0.37 ±0.09 


BSW I 


\H \ 2 




0.46 ± 0.04 


0.001 ±0.006 


0.43 ±0.05 


0.03 ±0.05 




\H-\ 2 




0.49 ±0.08 


0.91 ±0.10 


0.57 ±0.05 


0.97 ±0.06 




BR xlO" 


-3 


0.80 ±0.17 


0.52 ±0.12 


0.71 ±0.16 


0.52 ±0.11 


BSW II 


\H \ 2 




0.39 ±0.06 


0.07 ±0.04 


0.36 ±0.06 


0.11 ±0.06 




\H-\ 2 




0.58 ±0.08 


0.89 ±0.07 


0.63 ±0.06 


0.87 ±0.10 




BR xlO" 


-3 


0.79 ±0.19 


0.43 ±0.11 


0.62 ±0.17 


0.36 ±0.09 


CDDFGN 


\H \ 2 




0.46 ±0.04 


0.001 ±0.006 


0.43 ±0.05 


0.031 ±0.05 




\H_\ 2 




0.49 ±0.08 


0.91 ±0.10 


0.57 ±0.05 


0.97 ±0.06 




BR xlO" 


-3 


0.79 ±0.19 


0.43 ±0.11 


0.62 ±0.17 


0.36 ±0.09 


AW 


\H \ 2 




0.46 ± 0.04 


0.001 ±0.006 


0.43 ±0.05 


0.031 ±0.05 




W-\ 2 




0.49 ±0.08 


0.91 ±0.10 


0.57 ±0.05 


0.97 ±0.06 




BR xlO" 


-3 


1.10 ±0.25 


0.56 ±0.15 


0.81 ±0.23 


0.47 ±0.12 


ISGW 


\H \ 2 




0.46 ± 0.04 


0.001 ±0.006 


0.44 ±0.05 


0.03 ±0.05 




\H.\ 2 




0.48 ±0.08 


0.90 ±0.10 


0.56 ±0.05 


0.97 ±0.05 


Experiment |16] 


B[B -> ip(2S)K*) = (0.9 ± 0.29) x lO" 3 



and both first and third solutions of the ISGW model) are consistent with available 
experimental results, (ii) The predictions of the BSW I, BSW II and CDDFGN 
models show very little sensitivity to the solution- type, (in) \H + \ 2 has a very small 
value in all the formfactor models considered. 



3.3 B -> i/j(2S)K* Decays 

The branching ratio, longitudinal polarization and transverse polarization (|ff_| 2 ) for 
B — > ip(2S)K* can also be related directly to those for B — > %pK* by eliminating £i, 
£21 and £vi- The result is shown in Table To the best of our knowledge, only the 
branching ratio for this process is available |jl6| , 



B(B -> i)(2S)K*) = (0.9 ± 0.29) x 10" 



(42) 



From Table ||] we notice that the predictions of the branching ratio, \Hq\ 2 and 
|i7_| 2 are almost model independent. The predictions of solutions 1 and 3 are the 
closest to available experimental data. The other two solutions yield low branching 
ratios and close to zero longitudinal polarization. 
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4 Results in Factorization Approximation 



For the sake of completeness, we present in this section the factorization approxima- 
tion (£1 = £21 = £vi = 1) prediction for branching ratio, \H \ 2 and \H-\ 2 using the five 
theoretical models considered for the formfactors. These predictions are presented in 
Table 0. 

From Table |7|, we see that the factorization approximation predicts low values 
for the branching ratios compared to experiment. For the processes B — ► ipK* and 
B s — > tp<f), with the exception of BS W I model, the factorization approximation also 
underestimates the longitudinal polarization. 

If we scale the branching ratio by a factor of 3.5, we find that the BS W I model 
predictions agree, within error, with the available experimental data. This could be 
achieved by giving the nonfactorization parameters the values, 

£ 1 = v ^5, £ 21 = £ yi = l (43) 

or 

£1 = 6 = = V^5 . (44) 
This is what has been referred to as new factorization in ||. 

5 Discussion 

We have shown that B — ► ipK and ip(2S)K data require nonfactorized contributions 
in all of the five formfactor models we have considered. The smallest amount of 
nonfactorized contribution is needed for the BSW II model while the AW model 
requires the largest. We have calculated the branching ratios for B s — > iprj, iprf , 
ip(2S)r], and ip(2S)r]' in each model as functions of £i^ and displayed the result in 
Fig. H These branching ratios averaged over the five models are tabulated in Table 
[5] for a few values of the parameter £ Fl . 

We have used the latest CLEO data on B(B —>■ ipK*) and the transversity am- 
plitudes to determine the three nonfactorization parameters £1, £2 and £y defined 
in ([19]). There are four solutions which are expressed in terms of £1 and the ratios 
£21 = £2/£i, and £yi = £v/£i- These solutions are displayed in Fig. || for BSW I 
model and in Table [| for all the five models. We find that solutions exist for £2,1 = 1 
and £y,i = 1 only in BSW I model but then £1 7^ 1, i.e. there are no solutions where 
£1 = £2 = £v = 1 which would signal factorization. 

Assuming that the parameters £21 and £yi determined from B — > ipK* are the 
same in B s — ► ipip decay, we calculated the branching ratio, longitudinal and trans- 
verse polarizations of B s — > ifxp decay in all five models. Present data are consistent 
with model predictions but for a few exceptions as is seen from Table |5|. 

The branching ratio, longitudinal and transverse polarizations of B — > ip(2S)K* 
decay were also calculated assuming that the parameters £21 and £yi determined from 
B — > ipK* are the same in B — > ip(2S)K* decay. The results are shown in Table |6|. 
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Table 7: Predictions of branching ratios, \Ho\ 2 and \H-\ 2 for the processes 
B — ► ^-fiT*, 5 S — > an d -B — > ip(2S)K* using the factorization approximation 
£1 — £21 — £vi — 1- The errors in branching ratios are due to the errors in Wilson 
coefficients, decay constants and B meson life times. 





BR x KT 3 


\H \ 2 


\H-\ 2 






BSW I 




0.40 ± 0.24 


0.57 


0.39 






BSW II 




0.33 ±0.20 


0.35 


0.64 


B - 




CDDFGN 




0.24 ±0.14 


0.37 


0.62 






AW 




0.33 ±0.19 


0.12 


0.87 






ISGW 




0.15 ±0.09 


0.06 


0.93 






Experiment |7] 




1.35 ±0.18 


0.52 ±0.08 


0.47 ±0.08 






BSW I 




0.27 ±0.16 


0.55 


0.41 






BSW II 




0.23 ±0.14 


0.35 


0.64 


B s - 




CDDFGN 




0.22 ±0.13 


0.32 


0.66 






AW 




0.44 ±0.26 


0.20 


0.79 






ISGW 




0.09 ±0.06 


0.20 


0.80 






Experiment |21, 22 1 


0.93 ±0.33 


0.56 ±0.21 








BSW I 




0.24 ±0.14 


0.49 


0.43 






BSW II 




0.23 ±0.13 


0.29 


0.69 


B - 




CDDFGN 




0.13 ±0.07 


0.35 


0.65 






AW 




0.22 ±0.13 


0.23 


0.76 






ISGW 




0.14 ±0.08 


0.22 


0.77 






Experiment 16| 




0.90 ± 0.29 
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Finally, we find that in the factorization approximation none of the formfactor 
models predict correctly the branching ratios for for B — > ipK*, B s — > ip(f) and 
B — > ip(2S)K*. As for the longitudinal polarizations, |-ffo| 2 , m B — > ipK* and 
B s — > ^0, only the BSW I model predicts them correctly. BSW I also predicts \H_\ 2 
correctly for B — > ipK*. 

This research was partly supported by a grant to A.N.K. from the Natural Sciences and 
Engineering Research Council of Canada. 
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